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Abstract 

In this paper we explore the possibility that (a) fermion masses, (b) neutrino oscillations, (c) CP 
non-conservation and (d) flavor violations get intimately linked to each other within supersymmetric 
grand unification based on SO(10) or an effective G(224) = SU(2) L x SU(2) R x SU(4) C symmetry. 
We extend the framework proposed previously by Babu, Pati and Wilczek (BPW) which successfully 
described fermion masses and neutrino oscillations, to include CP violation. Assuming flavor universal 
SUSY breaking parameters at M* Z Mgut, and that CP violation arises through phases in the fermion 
mass matrices, we show how the presence of GUT threshold induces new and calculable CP and flavor 
violations. Including SM and SUSY contributions, we find that the BPW framework can correctly 
account for the observed flavor and/or CP violations in ArriK, Ams d , S(Bd — ► J/ipKs) and £k- While 
SUSY-contribution is small few%) for the first three quantities, that to 6k is sizable (~ 20-25%) 
and negative (as desired) compared to that of the SM. The model predicts S(Bd — > <f}Ks) to be in the 
range +(0.65-0.73), close to the SM prediction. The model yields Re(e'/e)s USY ~ +(4 - 14) x 10~ 4 ; 
the relevance of this contribution can be assessed only when the associated matrix elements are known 
reliably. The model also predicts that the electric dipole moments of the neutron and the electron should 
be discovered with improvements in current limits by factors of 10 to 100. 
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1 Introduction 



Ongoing studies of CP and flavor violations may well turn out to provide some important 
clues to physics beyond the standard model (SM). Strong agreement with the expectations 
of the SM would provide severe constraints on new physics, while departures could suggest 
the nature and the threshold for such new physics. On the experimental side there are now 
four well measured quantities reflecting CP and/or AF = 2 flavor violations. They are: 1 



where S(B d — > J/^>K S ) denotes the asymmetry parameter in (B d versus B d ) — > J/^>K S 
decays. It is indeed remarkable that the observed values including the signs of all four 
quantities as well as the empirical lower limit on ArriB a can consistently be realized within 
the standard CKM model for a single choice of the Wolfenstein parameters [1]: 



This fit is obtained using the observed values of e K = 2.27xl0~ 3 , V us = 0.2240 ± 0.0036, 
\V ub \ = (3.30 ± 0.24) xlO- 3 , \V cb \ = (4.14 ± 0.07) xl0~ 2 , |Am B J = (3.3 ± 0.06) xl0~ 13 GeV 
and Am Bd / AmB s > 0.035, and allowing for uncertainties in the hadronic matrix elements 
of up to 15%. One can then predict the asymmetry parameter S(B d — > J/^K S ) in the 
SM to be ~ 0.685 ± 0.052 [1,2]. This agrees remarkably well with the observed value 
S(B d ->• J/^K s )ex P t. = 0.734 ± 0.054, representing an average of the BABAR and BELLE 
results [3]. This agreement of the SM prediction with the data in turn poses a challenge for 
physics beyond the SM, especially for supersymmetric grand unified (SUSY GUT) models, 
as these generically possess new sources of CP and flavor violations beyond those of the SM. 

The purpose of this paper is to address the issues of CP and flavor violations 
in conjunction with those of fermion masses and neutrino oscillations, in the context of 
SUSY grand unification. In particular, our goal would be to obtain a unified description, 
in accord with observations, of all four phenomena: (i) CP non-conservation, (ii) flavor vio- 
lation, (iii) masses and mixings of quarks and leptons, as well as (iv) neutrino oscillations, 
within a single predictive framework based on SUSY grand unification. 



x t' K reflecting direct AF = 1 CP violation is well measured, but its theoretical implications are at present unclear due to 
uncertainties in the matrix clement. We discuss this later. 



Am K , e K , Am Bi and S(B d -> J/^K S ) 



(1) 



p w = 0.178 ± 0.046; fj w = 0.341 ± 0.028 . 



(2) 



2 



In this regard, (a) the observed quantum numbers of the members of a family, (b) gauge 
coupling unification, (c) neutrino oscillations, as well as (d) the likely need for leptogenesis [4] 
as a prelude to baryogenesis, together, suggest that the SM gauge symmetry very likely 
emerges, near the GUT scale Mjj ~ 2 x 10 16 GeV, from the spontaneous breaking of a 
higher gauge symmetry that possesses the symmetry SU(4)-color [5]. The higher gauge 
symmetry in 4D could be either SO(10) [6], or an effective (presumably string derived) 
symmetry G(224) = SU(2) L x SU(2) R x SU(A) C [5]. The need for SU(4)-color arises because 
it provides: (a) the right-handed neutrino (vr) as a compelling member of each family, 
(b) B — L as a local symmetry, and (c) the Dirac mass of v T in terms of the top-quark 
mass. These three ingredients, together with the SUSY unification scale, seem to be rather 
crucial [7] to an understanding of the neutrino masses via the seesaw mechanism [8] and 
implementing baryogenesis via leptogenesis [4] . The theory thus described may be viewed to 
have emerged in 4D from a string/M theory near the string scale M st > M GUT 2 . It should 
of course possess weak scale supersymmmetry so as to avoid unnatural fine tuning in the 
Higgs mass and ensure gauge coupling unification. 

A predictive framework based on the symmetry SO(10) or G(224), and a minimal Higgs 
system which accomplish all these has been proposed in Ref. [9], which we refer to as the 
BPW model. This model describes the masses and mixings of all fermions including neutrinos 
by making the simplifying assumption that the fermion mass matrices are real and thus 
CP-conserving. Notwithstanding this assumption, the framework is found to be remarkably 
successful. In particular, it makes seven predictions involving fermion masses, CKM elements 
and neutrino oscillations, all in good accord with observations, to within 10% (see discussion 
in Sec. 2). 

Now in general one would of course expect the entries in the fermion mass matrices 
to have phases because the VEVs of the relevant Higgs fields, and/or the effective Yukawa 
couplings, can well be complex. These in turn can induce CP violation through the SM/ CKM 
interactions as well as through SUSY interactions involving sfermion/gaugino loops. 

The question arises: Can the BPW- framework of Ref. [9], based on the supersymmet- 
ric SO(10) or G(224)-symmetry, be extended, by allowing for phases in the fermion mass 

2 The relative advantages of SO(10) versus G(224) as effective 4D symmetries involving the issues of coupling unification on 
the one hand and doublet-triplet splitting on the other hand are discussed in [7]. 
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matrices, so as to yield net CP and flavor-violations, arising through both SM and SUSY in- 
teractions, in accord with observations, while still preserving its successes as regards fermion 
masses and neutrino oscillations? This is the question we pursue in this paper. 

As we will see, these four phenomena - (i) fermion masses, (ii) neutrino oscillations, (iii) 
CP non-conservation, and (iv) flavor violations - get intimately linked to each other within 
the SUSY SO(10)/G(224) framework. Satisfying simultaneously the observed features of all 
four phenomena within such a predictive framework turns out, however, to be a non-trivial 
challenge to meet. One main purpose of the present paper is to show that the answer to 
the question raised above is in the affirmative. We defer the discussion of our results and 
predictions in this regard to section 5. 

Since CP violation would have its origin, in our model, entirely in the fermion mass 
matrices, and since the BPW framework presented in Ref. [9] has proven to be so successful 
in describing fermion masses and neutrino oscillations, we first briefly recall the salient 
features of this prior work in the next section. In the following section we extend the same 
so as to include CP violation. 

2 Fermion Masses and Neutrino Oscillations in G(224)/SO(10): 
A Brief Review of Prior Work 

The 3x3 Dirac mass matrices for the four sectors (u, d,l,u) proposed in Ref. [9] were 
motivated in part by the notion that flavor symmetries [10] are responsible for the hierarchy 
among the elements of these matrices (i.e., for "33">"23">"22">"12">"11", etc.), and 
in part by the group theory of SO(10)/G(224), relevant to a minimal Higgs system. Up to 
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minor variants [11], they are as follows: 



M D = 





-e' 




3e' 
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These matrices are defined in the gauge basis and are multiplied by ty L on left and 
on right. For instance, the row and column indices of M u are given by (ml,^,^) and 
(ur, Cr, £r) respectively Note the group-theoretic up-down and quark-lepton correlations: 
the same a occurs in M u and MjP, and the same rj occurs in Md and Mi. It will become clear 
that the e and e' entries are proportional to B — L and are antisymmetric in the family space 
(as shown above). Thus, the same e and e' occur in both (M u and M d ) and also in (MjP and 
Mi), but e — > — 3e and e' — > — 3e' as g — > /. Such correlations result in enormous reduction 
of parameters and thus in increased predictiveness. Such a pattern for the mass-matrices 
can be obtained, using a minimal Higgs system 45#, 16#, 16// and 10# and a singlet S of 
SO (10), through effective couplings as follows [12] (see Ref. [9] and [7] for details): 

Auk = ^3316316310^ + [/i 2 316 2 16 3 10h(S , /M) 

+ a 23 16 2 16 3 10 H (45 H /M')(S/MY + ( 723 16 2 16 3 16^(16h/M /, )( 1 S/M)' ? ] 

+ [h 22 1 6 2 1 6 2 1 H (S/M) 2 + ( ?22 1 6 2 1 6 2 1 6^(16 ff /M /, )( 1 S/M)' ?+1 ] 

+ [ 9l2 lQ 1 lG 2 l§ d H {lQ H /M"){S/My +2 + a 12 16 1 16 2 10 H (45 H /M')(S/M) p+2 ] . (4) 

Typically we expect M', M" and M to be of order M str i ng or (possibly) of order Mqut [13]. 
The VEV's of (45//) (along B — L), (16//) = (16//) (along standard model singlet sneutrino- 
like component) and of the SO(10)-singlet (S) are of the GUT-scale, while those of 10// and 
of the down type SU(2)/,-doublet component in 16// (denoted by 16#) are of the electroweak 
scale [9,14]. Depending upon whether M'(M") ~ M GUT or M string (see [13]), the exponent 
p(q) is either one or zero [15]. 

The entries 1 and a arise respectively from /i 33 and h 23 couplings, while fj = rj — a and rf 
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arise respectively from g 2 3 and (^-couplings. The (B — L)-dependent antisymmetric entries 
e and e' arise respectively from the 023 and a i2 couplings. [Effectively, with (45 h) oc B — L, 
the product 10# x 45h contributes as a 120, whose coupling is family-antisymmetric] The 
small entry arises from the ^22-coupling, while C22 arises from the joint contributions of 
h 2 2 and ^22-couplings. As discussed in [9], using some of the observed masses as inputs, one 
obtains \fj\ ~ \a\ ~ |e| ~ (9(1/10), \r]'\ « few x 10~ 3 and |e'| ~2x 10~ 4 . The success of the 
framework presented in Ref. [9] (which sets £22 = C22 = 0) in describing fermion masses and 
mixings remains essentially unaltered if | (C225 C22) I — (1/3)(10 -2 ) (say). 

Such a hierarchical form of the mass-matrices, with /t 33 -term being dominant, is attributed 
in part to flavor gauge symmetry(ies) that distinguishes between the three families, and in 
part to higher dimensional operators involving for example (45n)/M f or (16#) /M" , which 
are suppressed by MGUT/M*rmg ~ 1/10, if M' and/or M" ~ M str i n g. As an example, 
introduce just one U(l)— flavor gauge symmetry, together with a discrete Z 2 — symmetry D, 
with one singlet S. The hierarchical form of the Yukawa couplings, in accord with Eqs. (jHj) 
and (HJ), would follow, for the case of p = 1, q = 0, if, for example, the U(l) flavor charges 
are assigned as follows: 3 

16 3 16 2 16i 10^ 16 H 16 H 45 H S 

(5) 

a a + 1 a + 2 -2a -a -1/2 -a -1 

All the fields are assumed to be even under the discrete symmetry D, except for 16# and 16 h 
which are odd. It is assumed that other fields are present that would make the U(l) symmetry 
anomaly-free. With this assignment of charges, one would expect |C22 d | ~ {(S)/M) 2 ] one 
may thus take |C 2 2 d | ~ (1/3) x 10~ 2 without upsetting the success of Ref. [9]. In the same 
spirit, one would expect |Cxa, Csx I ~ ((S)/Mf ~ 10" 2 and |Cu| ~ ((S)/M) 4 ~ 10" 4 (say), 
where (n, C13, and C31 denote the "11", "13", and "31", elements respectively. The value of 
a can get fixed by the presence of other operators (see Ref. [16]). 

The Majorana masses of the RH neutrinos arise from effective couplings of the form 
Aviaj = /ijl6il6jl6i/16#/jW, where the /y's include appropriate powers of (S/M) in accord 

3 An alternative assignment with p = 1, q = 1 would arise by choosing the charges: (0,1,2,0,0,0,0 and —1) for 
(I63, 162 , 161 , IOjj, 16/j, 16h , 45h and S) respectively. Another variant would be to assign a flavor charge of 3 (or a+3 
in case of Eq. JSJ) to the first family I61, leaving all other charges as above. This would suppress "13", "31", "12", "21", and 
"11" entries by an order of magnitude relative to the other cases. 
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with flavor-charge assignments of 16j and 16h, and M is expected to be of order string or 
reduced Planck-scale. From the flavor-charge assignments given in Eq. (J5|). it is clear that 
the "33" entry is leading and the other entries in the Majorana mass matrix have a relative 
hierarchy which is identical to that in the Dirac mass matrices. We refer the reader to [9] 
and [7] for a more detailed discussion of the neutrino sector. Here, we list only the results. 

Ignoring possible phases in the parameters and thus the source of CP violation for 
a moment, and also setting = C22 = O5 as was done in Ref. [9], the parameters 
(a, 77, e, e', 7/, Ai^, and can be determined by using, for example, m\ hys = 174 GeV, 
m c (m c ) = 1.37 GeV, m 8 (l GeV) = 110 - 116 MeV, m u (l GeV) = 6 MeV, and the observed 
masses of e, /1, and r as inputs. One is thus led, for this CP conserving case, to the following 
fit for the parameters, and the associated predictions [9]: 

a « 0.110, 77 « 0.151, esi -0.095, \r]'\ « 4.4 x 10" 3 , 

(6) 

e' « 2 x 10- 4 , Ml « m t (M x ) « 100 GeV, M° d « m r (M x ) « 1.1 GeV. 

These output parameters remain stable to within 10% corresponding to small variations 
(~ 10%) in the input parameters of m t , m c , m s , and m u . These in turn lead to the following 
predictions for the quarks and light neutrinos [9], [7]: 



m b {m b ) « 4.9 GeV, y/Arn^ « m(u 3 ) w (1/24 eV)(l/2-2), V cb w 0.044, 
sin 2 2^ T w 0.98 - 0.995, (for m(v 2 )/m{v z ) w 1/10 - 1/5), (7) 
IK. J w 0.20, 



0.07, m d (l GeV) w 8 MeV 

It has been noted [7, 16] that small non-seesaw contribution to v\y^ mass term (~ few 
xlO -3 eV) which can arise through higher dimensional operators, but which have been 
ignored in the analysis given above, can lead quite plausibly to large v e — oscillation 
angle in accord with the LMA MSW solution for the solar neutrino problem. Leaving aside 
therefore the question of the v e — oscillation angle, it seems quite remarkable that all 
seven predictions in Eq. (JJJ) agree with observations to within 10%. Particularly intriguing 
is the (B — L)-dependent group-theoretic correlation between V cb and 0° sc Vt , which explains 
simultaneously why one is small (V cb ) and the other is so large (9° S 1 T ) [7,9]. That in turn 
provides some confidence in the gross pattern of the Dirac mass matrices presented above 
and motivates the study of CP and flavor violations within the same framework. This is 
what we do in the next section. 
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3 Phases in the Fermion Mass Matrices: The Origin of CP viola- 
tion 



In the work of Ref. [9] reviewed above, the parameters (a, e, 77, e', r( etc.) entering into 
the fermion mass matrices were assumed to be real, for simplicity, and thereby (at least) the 
SM interactions were rendered CP-conserving 4 . Noting that the VEVs of the Higgs fields 5 
and/or the effective Yukawa couplings can well be complex, however, we now propose to 
extend the SO(10)/G(224) framework reviewed above to include CP violation by allowing 
for these parameters to have phases. 

Given the empirical constraints on (i) CP and flavor violations, as well as (ii) fermion 
masses and (iii) neutrino oscillations, on the one hand, and (iv) the group-theoretical con- 
straints of the SO(10) / G(224) framework on the other, it is of course not at all clear, a priori, 
whether any choice of phases and variations in the parameters of the fermion mass matrices 
presented above can yield observed CP and flavor-violations, and simultaneously preserve 
the successes of the framework of [9] as regards fermion masses and neutrino oscillations. 
That is the issue we now explore. We choose to diagonalize the quark mass matrices M u and 
Ma at the GUT scale ~ 2 x 10 16 GeV, by bi-unitary transformations - i.e. 

M diag = X d\ MdX d R M di ag = X f M U X% (8) 

with phases of q l L R chosen such that the eigenvalues are real and positive and that the CKM 
matrix V C km (defined below) has the Wolfenstein form [17]). Utilizing the hierarchical nature 
of the mass matrices, one can obtain (approximate) analytic expressions for the diagonalizing 
matrices. They are: 





e -i((p v - e ) 


Irf/Xdle-tfo-'+M 


fj'^fj — g|g *(<^»7-<! C33) 


xi~ 


-\q' / X d \e i{f h + < +4 ' x J 


e i(4> v+e +4>x d -C,us) 


_|_ g | e *(0T)+e+C33) 




\r]' / X d \\r] + e^xJ -Y 


-\r] + e\e i((l)x <i-^ 





4 modulo the contribution from the strong CP parameter 



5 For instance, consider the superpotential for 45/i only: W(45h)= M4s(45j/) 2 + "A(45/-f) 4 "/M, which yields (setting 
^45^=0), cither (45j-j) = 0, or (45^) 2 = — (2M45M/ "\"). Assuming that "other physics" would favor (45jy) ^ 0, we see that 
(45 jf) would be pure imaginary, if the quantity in the brackets is positive with all parameters being real. In a coupled system, 
it is conceivable that (45 h) m turn would induce phases (other than and 7r) in some of the other VEVs as well, and may 
itself become complex rather than pure imaginary. 



8 



X 



R 



'|r] / /Xd|e _i< - 9!>, '" £ - ) g-i(^7,-e+C«s) |^ _ e | e — *C<^tj— e— C33) 



|r/7X d ||r7 - e l 



— \r] — e\e 



-iC,u 



e*?33 



(10) 



Here 4> v ±e = arg(rj ± e), that is, (77 ± e = [77 =J= e| e l( ^> ±E ); K = 77' [ 77 — e|e ^ ud and 



X d = -|77 2 - e 2 | + |&| e " 



Pv+€+<Pv-^-'P(2d 



] = \X d \e i,f>x d, where C22 = IC2 2 l e ^ C2<i - The corre - 



sponding matrices X"l R for diagonalizing the up sector are obtained from above with the 
substitutions : 77 — > er; (22 ~ > C22! (v'^ 6 ') ~^ ^ ze '- Thus (j) n ± e are replaced by <f) a ± e = arg(cr±e)] 
and X d by X u = -\a 2 - e 2 | + |^|e- < to'+« + ^-« - *feJ = |X u |e ife «. Given the definitions of 



k >x d and cj)x u as above, we have 



=33 



'ij+e, 



+ R sin Q; 7 = (0^ + </>,,_,) - + </> CT _ e ) + 



where 



R = \e'/X u \/\r]'/X d \ ps \/m u /m c ^/m s /m d « 0.3; 

/3 n « i2(sinfi/n), n = (0x d - 0xj + 7; 

C cb ~ ar^fe^^d^ + el- I t7 + e|e <( ^ + «"^ + ' ) }]; 

Cms ~ — -R sinf2[l — i? cosfi] -1 . 



As mentioned above, using observed fermion masses and mixings [9], we obtain: |e'| ~ 
1/10 |?/|, with \rj'\ ~ (few)xl0~ 3 < (\rj\ ~ |e| ~ \a\ ~ 1/10) < I. In writing Eqns. © 
and (jlOj) . we have not displayed for simplicity of writing, small correction terms {0(e 2 , 1] 2 )), 
which are needed to preserve unitarity. We have also not displayed small phases of order 
\r]'e'/X u Xd\ X sin Q ~ 1/100, \e'/r]'\ ~ 1/10 and RsinQ ~ 1/10. Our results to be presented, 
that are based on exact numerical calculations, however incorporate these small corrections. 

The CKM elements in the Wolfenstein basis are given by the matrix Vckm — 



e- ia {Xl ] Xl), where a = {<f)^ t 
(given 1 7/ I |e'|), we can choose </y. 



■"rj—e ) 



- (0 e / — <f) v i +e r), where without loss of generality 
0. To a good approximation, the CKM elements 
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are given by: 

V ud » V cs *V tb nl 

V us »| |7//X d | - |e7X u |e in | » 

V cb «| e^-^^l + e | - | cx + e | e^' - *^} | « -V^ 

Kb ~ [V | ?7 - e I - I e'/X u I e^-^{|?7 + e| - |<r + e | e -^+^-^)}] x e^a+AiJ-C*] 
Vid « [| r^'/Xrf I e i(fe d){| e + r] I - I o- + e I e -*(^+«-*»+«)} _ j/^ _ e |] x e -*[«(i+fti)-Cc6] 

Note that the CKM elements have the desired Wolfenstein form with only V ub and Vtd being 
complex and the others being real to a good approximation. defined above is just the 
argument of the expression within the bars for V cb . One can check that to a good approxima- 
tion, (neglecting the r]'\r]—e\ term for V td that causes < 10% error), the phase of V td is given by 
<p td = ArgiVu) ss -RsmVt, and (V^l w IV/^dllK*! ~ \Jm d /m s \V ch \, and similarly 
|Kb| ~ \/m u /m c \V cb \. 

Before presenting the results of a certain fit and the corresponding predictions, we need 
to first discuss SUSY CP and flavor violations in the presence of phases in the fermion mass 
matrices. This is done in the next section. 

4 SUSY CP and Flavor Violations 

Our procedure for dealing with SUSY CP and flavor violations may be summarized by the 
following set of considerations: 

1) As is well known, since the model is supersymmetric, non-standard CP and flavor 
violations would generically arise in the model through sfermion/gaugino quantum loops 
involving scalar {mass) 2 transitions. The latter can either preserve chirality (as in q l L R — > 
q 3 L R ) or flip chirality (as in q l L R — > q RL ). Subject to our assumption on SUSY breaking 
(specified below), it would turn out that these scalar {mass) 2 parameters get completely 
determined within our model by the fermion mass-matrices, and the few parameters of 
SUSY breaking. 

2) SUSY Breaking : We assume that SUSY breaking is communicated to the SM sector 
by messenger fields which have large masses of order M*, where M GUT < M* ps M string , such 
that the soft parameters are flavor-blind, and family- universal at the scale M*. A number 
of well motivated models of SUSY breaking, e.g., those based on mSUGRA [18], gaugino- 
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mediation [19], anomalous U(l) — D term [20,21], combined with dilaton-mediation [21,22], 
or possibly a combination of some of these mechanisms, do in fact induce such a breaking. 
While for the first two cases [18,19] we would expect extreme squark degeneracy (ESD) i.e. 
k = \m 2 (qi) — rn 2 (qj)\/m 2 (q)Av ^ 10~ 3 (say) at the scale M*, for the third case [20,21], 
one would expect intermediate squark degeneracy (ISD) i.e. k ~ 10~ 2 (1 — 1/3) at M*. 
For the sake of generality, we would initially allow both possibilities, k — (ESD), and 
k ~ 10- 2 (1 - 1/3) (ISD) at M*. 

In an extreme version of universality, analogous to CMSSM, the SUSY sector of the model 
would introduce only five parameters at the scale M*\ 

m o, m i/2, A Q , tan/5 and sgn(n). 

In some cases, A Q can be zero or extremely small (< lGeV) at M* as in [19] and [21]. 
For most purposes we will adopt this restricted version of SUSY breaking, including the 
vanishing of A at M*. However, our results will be essentially unaffected even if A is non- 
zero (~ 500 GeV, say) but real (see remarks later). We will not insist on, but will allow for, 
Higgs-squark-slepton universality, which does not hold, for example, in the string-derived 
model of [21]. In spite of flavor-preservation at a high scale M* , SUSY-induced flavor- 
violation would still arise at the electroweak scale through renormalization group running of 
the sfermion masses and the A-parameters from M* — > Mqut ^w, as specified below. 
Although the premises of our model as regards the choice of universal SUSY parameters 
coincide with that of CMSSM, as we will see, owing to the presence of GUT-scale physics 
in the interval M* — > Mqut, SUSY CP and flavor violations in our model (evaluated at the 
electroweak scale) would be significantly enhanced compared to that in CMSSM (or even 
CMSSM with right-handed neutrinos). This difference provides some distinguishing features 
of our model. 

3) Flavor Violation due to RG Running of Scalar Masses from M* to M GU t 

For MSSM embedded into SO (10) above the GUT scale, there necessarily exist heavy 
color-triplet Higgs fields which couple to fermions through the coupling h t 16316310^, while 
there exist heavy doublets for both SO(10) and G(224) which also couple to fermions owing 
to the mixing of 10 h with 16 h (see [9]). (Here h t stands for /i 33 of Eq. (@J)). These couple 
to and 6r with the large top quark Yukawa coupling h t . The heavy triplets and doublets 



11 



possess masses of order M GUT . One can verify (see [23]) that the evolution of RG equations 
for squark masses involving such couplings suppress bh and Sr masses significantly compared 
to those of d^R and §l,,r- Note that left-right symmetry implies equal shifts in bi and 6r 
masses arising from GUT scale physics in the momentum range Mqut < < M*. Such 
differential mass shifts i.e.- (mj — rh\ ^ L R = Am- , for the embedding of MSSM into 
SO(10), are found to be (with A Q = 0): 

30m 2 

Ami = A ™i « -(T#)^ HM*/M GUT ) = -(ml/m . (12) 

The hat signifies GUT scale values. Here m Q denotes the (approximately) degenerate mass of 
squarks at the scale M*. We have set hf = 1/2; we expect M*/Mqut ~ (3 to 10), say, and 
thus, f = ln(M*/M GUT )/2.6 w (0.4 to 0.9). For the case of MSSM embedded into G(224), 
which provides the heavy doublet, but not the triplets, the factor 30 in Eq. (|12[) should be 
replaced by 12. 

Having diagonalized the quark mass-matrices and Mu^ at the GUT scale by matrices 
as in Eq. (JHJ), SUSY flavor violation may be assessed by imposing the parallel transformations 
on the squark (mass) 2 matrices {(M^) ll/rr) defined in the gauge basis, i.e., by evaluating 
(M^)llX£ and X^(M^)rhX r , and similarly for the up sector. Following discussion 
on SUSY breaking, the off-diagonal elements (in the gauge basis) and all chirality flipping 
elements are set to be zero - i.e {M^^ll/rr = (i ^ j) and (A^)lr = - at the scale M*. 
Once squarks are non degenerate at M G ut owing to the mass-shift of bi,R as in Eq. (|T2*|) . 
the transformations mentioned above induce off-diagonal elements with squarks being in the 
SUSY basis. For the down squark mass matrices (evaluated at the GUT scale), these off 
diagonal elements are found to be: 
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8f L {M GUT ) ~ [Kf SD + (Am? £ /my {-\i/X d \\e + r/| 2 + ?/|e 2 - r/ 2 |e^)}] e -^ 

<&(M GC/T ) ^ [k? 3d + (Aml R /m%)(-\ V '/X d \\e - V \ 2 + V '\e 2 - r, 2 \e~^)}} e -^ 
« [k}| d + 3£ x 1(T 3 - 10- 5 (Oe"^](m 2 /m 2 9 )e-^ d 

<5f L (M Gl/T ) ~ (Amf /m 2 g ) [-jj'ltj - e|e*3s + \r)'/X d \\r) + ele^V^)] 
« [(2.50 x 10- 4 e^3 - (2.50 x l^e^-^iml/m 2 ^) 

(13) 

S R 3 R (M GUT ) ~ (Am^/m 2 g )[-V|r? + e^s-^) + |,//JSQ| |?7 - e|e iC 3 3 ] 
« -[(1.250 x 10^ 2 e i? 3 d 3 ]( m 2/ m 2j 



« [(1.250 x 10- 2 e i(c 33-^+^)]( m 2/ m 2j 



- (Afhf /my - e\e^+M 



w [(6.20 x 10~V( c 33+<^)]( m 2 /m 2 

Here ^ 2 D = [(mS° )2 - m? )2 )]/m 2 g (|VAd|) ~ ±(2 x 1(T 3 )(1 - 1/3); this term would 
be present for the case of intermediate squark degeneracy (ISD), corresponding to small 
(~ 10~ 2 (1 — 1/3)) squark non-degeneracy at the scale M*, as in models of Ref. [20,21]. From 
now on, for the sake of concreteness, we drop this term, 6 setting k 12 sd = 0. In above <ptd ~ 

- \e'/X u \/\r]'/X d \ sinfi ~ (-1/3) sinfi ~ (-1/6) (say). The hat on top signifies GUT 
scale values, and 8l L / rr = (A^ i / iJK )/m 2 g , where K l [ L denotes the (mass) 2 parameter for 
QjL Qn transition in the SUSY basis. Here, m sq denotes the average mass of the di,,R 
and sl,r squarks, which remain nearly degenerate( to 1% or better) even at the weak scale. 
For each 8 l l L i RR we have exhibited approximate numerical values by inserting values of the 
parameters r], a, e rf etc. for some typical fits (as in Eqs. © and (|2*0jl) to indicate their 
typical values. 

"Note that the case of ISD {k) 2 sd ~ (2 X 10~ 3 )(1 - 1/3) j= 0) would make a difference only for the case of K° — K° transitions 

- that is, for Amj; and ejf. 
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Assuming for simplicity, universality of scalar masses m Q (of the first two families) and of 
the gaugino masses m.1/2 at the GUT scale, the physical masses of squarks of the first two 
families and of the gluino are given by: 

ml q & m 2 + 7.2 m 2 1/2 ; « 2.98 m 1/2 . (14) 

This result is rather insensitive to the mass shifts of bi,R- Using the above relations we get 
p x = (m 2 Jm 2 sq ) ~ 1 - 0.8a; » (0.84, 0.76, 0.5 and 0.2) for x = mj/m 2 sq =(0.2, 0.3, 0.6 and 1), 
which enters into all the (^-elements in Eq. (II 3j) . 

We remind the reader that the elements 5 % [ L RR , induced solely through GUT scale physics 
being relevant in the interval M* — > Mgut, would be absent in a general CMSSM or MSSM, 
and so would the associated CP and flavor violations. 

4) Flavor Violation Through RG Running From M G ut to m w in MSSM : It is 
well known that, even with universal masses at the GUT scale, RG running from Mqut 
to mw in MSSM, involving contribution from the top Yukawa coupling, gives a significant 
correction to the mass of b' L = Vtddh + V ts si + VtbbL, which is not shared by the mass-shifts 
of bji, (Il^r and Sl,r- This in turn induces flavor violation. Here, d^, sl and bi are the SUSY 
partners of the physical di, sl and bi respectively. The differential mass shift of b L arising 
as above, may be expressed by an effective Lagrangian [24]: AC = — (Am' 2 )b£b L , where 7 

Am' 2 = -3/2m 2 T] t + 2.3A m 1/2 r] t (l - Vt ) - (A 2 /2) Vt (l - Vt ) + m 2 1/2 (3r] 2 - 7 Vt ) . (15) 

Here rj t = (h t /hf) w {m t /v sin/?) 2 (l/1.21) w 0.836 for tan/3 = 3. Numerically, setting 8 
A = 0, Eq. (JT3J) yields: (Am' 2 /m 2 q ) w -(0.40, 0.34, 0.26, 0.22) for x = mj/m 2 sq w 
(0.1, 0.4, 0.8, 1.0). Expressing b' L in terms of down-flavor squarks in the SUSY basis as 
above, Eq. (|T5|l yields new contributions to off diagonal squark mixing. Normalizing to m 2 q , 
they are given by 

C' 13 ' 23) = (^r)(m, Kv», v*v tb ) . (16) 

The net squark (mass) 2 off-diagonal elements at mw are then obtained by adding the 

7 Note that strictly speaking Eq. 1151 holds if the soft parameters are universal at the GUT-scale. However, the correction 
to this expression due to RG running from M* to mw would be rather small, being a correction to a correction. 

8 Although we have put A = (for concreteness), note that Arn^ 2 would typically get only a small correction (2; 5%), even 
if A were non-zero (< 1 TeV), with m sa (0.7 - 1) TeV and m 1/2 ss (200 - 300) GeV, say. 
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respective GUT-scale contributions from Eqs. (|13j) to that from Eq. They are: 

= & + $r = %r (17) 

From the expressions given above (Eqs. (fl"3"j) and (fT6*jO . it follows that for a given choice 
of the SUSY-parameters (i.e. m Q , m^j or equivalently m sq and rrig), SUSY CP and flavor 
violations are completely determined within our model by parameters of the fermion mass- 
matrices. This is the reason why within a quark-lepton unified theory as ours, SUSY CP 
and flavor violations get intimately related to fermion masses and neutrino oscillations. 

5) A — Terms Induced Through RG Running from M* to Mgut: Even if A D is 
zero at M* (as we assume, for concreteness, see also [19] and [21]), RG running from M* 
to Mqut i n the context of SO(10)/G(224) would still induce non-zero A parameters at the 
GUT scale [23]. For our case, the A terms are induced through loop diagrams involving 
the /133, #23, and 023 couplings and the SO(10) or G(224) gauginos. We find that if we take 
^10^ ~ Mw H ~ Mqut, we can write the Aj^-matrix at the GUT-scale for the down 
squark sector in the SUSY basis for the case of SO(10) as follows: 



A a LR = Zx (XI 



95e' + 90r/ 

-95e' + 90r/ 63(2*2 95e + 90r/ - 27a 

-95e + 90r] - 27a 63 



where Z = ( T ^)/i t (7 1 2 M A /n(-^), and (X d ) L)R are given in Eqs. © and (HJJJ) . The 323 
coupling does not contribute to the up-sector; thus the A-matrix for the up squarks, A^ R , 
can be obtained from above by setting 7/ = and replacing 90?7 — 27a by 63<r, (22 by 
(22; an d Xf R by XI R in A d . Similarly, the lepton ^-matrix, A l LR is obtained by letting 
(e, e') — > — 3(e, e') and replacing X d R by X LR in A LR . For the case of G(224), the matrix 
A d LR would be obtained by making the substitutions: (90, 63, 95) -> (42, 27, 43) in Eq. (fTHl) . 
and likewise in A U LR and A l LR . It is sometimes convenient to define the sfermion transition 
mixing angles as 

(Sfnh - (At l R h (^r) = (At l R h (Stnh - (A u LR h (£j) • (19) 

Note that these induced A^-terms for all three sectors, like the squarks (mass) 2 elements 
^llrr gi ven i n Eqs. (|13|) - (jl7j) . are completely determined within our model by the fermion 
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mass matrices, for a given choice of M\ w mi/2 and ln(M* / 'Mqut)- We now utilize these 
SUSY CP and flavor-violating elements to predict the results of our model. 

Once again, as in the case of o[ L RR , these induced A-terms arising purely through GUT- 
physics, would be absent or negligibly small in CMSSM. As a result, some of the interesting 
predictions of our model as regards e' K and edm's (to be discussed below) and lepton fla- 
vor violations (to be discussed in a forthcoming paper [25]) would be absent altogether in 
CMSSM. 

5 Compatibility of CP and Flavor Violations with Fermion Masses 
and Neutrino Oscillations in SO(10)/G(224): Our Results 

It has been noted in Sec. 1 that (given about 15% uncertainty in the matrix elements) the 
SM agrees very well with all four entries of Eq. (JT}, for a single choice of the Wolfenstein 
parameters pw and fjw (Eq. ©)• The question then arises (as noted in Sec. 1): If a SUSY 
SO(10) or G(224) model is constrained by requiring that it should successfully describe 
fermion masses and neutrino oscillations (as in Sec. 2), can it still yield (for some choice of 
phases in the parameters r], a, e etc.) values for p w and fjw more or less in accord with the 
SM-based phenomenological values for the same, as listed in Eq. (j2J? Anticipating that (for 
any given choice of the parameters rj, a, e etc.) the SO(10)/G(224) model-based values of 
Pw and fjw would generically differ from the SM-based phenomenological values (given in 
Eq.©), we will denote the former by p' w and fj' w and the corresponding contributions from 
the SM-interactions (based on p' w and fj' w ) by SM' . The question that faces us then is this: 
When the SM' contributions are added to the SUSY contributions arising from the three 
sources listed in Sec. 4, can such a constrained SO(10) or G(224) model account for the 
observed values of all the four quantities listed in Eq. and in addition is it consistent 
with the empirical upper limits on the edm's of the neutron and the electron? 9 . 

Before presenting our results, we make some preliminary remarks. First of all one might 
have thought, given the freedom in the choice of phases in the parameters of the mass 
matrices, that it ought to be possible to get almost any set of values of {pw and fjw), and 
in particular those in accord with the SM values (Eq. Q). It turns out, however, that in 

9 We extend the same question to include lepton flavor violating processes (such as /i — > e 7 and r — > [i 7) in a separate note. 
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general this is indeed not possible without running into a conflict with the fermion masses 
and/or neutrino oscillation parameters within a SO(10) or G(224)-model 10 . In other words, 
any predictive SO(10) or G(224)-model is rather constrained in this regard. 

Second, one might think that even if the SO(10)/G(224) model-derived entities p' w and 
fj' w , constrained by the pattern of fermion masses and neutrino oscillations, are found to be 
very different in signs and/or magnitudes from the SM values shown in Eq. (J2J), perhaps the 
SUSY contributions added to the SM' contributions(based on p' w and fj' w ) could possibly 
account for all four quantities listed in Eq. It seems to us, however, that this is simply 
not a viable and natural possibility, unless one is willing to invoke MSSM and finely adjust its 
arbitrary (in general some 105) parameters, as needed. In the latter case, the good agreement 
between experiments and the SM predictions would appear to be fortuitous (see Sec. 1). 

This is why it seems to us that the only viable and natural solution for any SUSY G(224) 
or SO(10) model for fermion masses and neutrino oscillations is that the model, allowing 
for phases in the fermion mass matrices, should not only yield the masses and mixings of 
all fermions including neutrinos in accord with observations (as in Sec. 2), but it should 
yield p' w and fj' w that are close to the SM values shown in Eq. (J2J). This, if achievable, 
would be a major step in the right direction. One then needs to ask: how does the combined 
(SM' + SUSY) contributions fare for such a solution as regards its predictions for the four 
quantities of Eq. JIJ and other CP and/or flavor violating processes, for any given choice of 
the SUSY parameters (m Q , mi/ 2 , A Q , tan/5 and sgn(p))! It should be stressed here that 
even if the CKM elements including p' w and ff w should turn out to be close to the SM values 
(Eq. (J2J), the SUSY contributions can in general still have a marked effect, in accord with 
observations, at least on some of the processes where the SM (or SM' )- contributions are 
naturally suppressed (as in the case for e#, edm's and lepton flavor violating transitions). 
Study of these processes, some of which we discuss below, can help distinguish between the 
SM versus the SUSY SO(10)/G(224)-models. 

Without further elaboration, we now present our main results. In this paper we will 
present only one fit to the parameters which has the desired properties. 11 Allowing for phases 
(~ 1/10 to 1/2) in the parameters 77, a, e', and ($2 °f the G (224) /SO (10) framework (see 

10 for a discussion of difficulties in this regard within a recently proposed SO(10)model, see e.g. [26] 
We have verified that there actually exists a class of fits which nearly serve the same purpose. Only one of these (Eq. 

is exhibited here for the sake of concreteness. 
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Eq. (J5J)) we find that there do exist solutions which yield masses and mixings of quarks and 
leptons, in accord with observations to within 10% for most part (see discussion below), and 
at the same time yield p' w and fj' w close to the SM values, as given in Eq. (J2J). A desired fit 
to the parameters is given by: 

a = 0.109 - 0.012i, n = 0.122 - 0.0464?;, e = -0.103, rf = 2.4 x 10" 3 , 

(20) 

e' = 2.35 x 10-V( 69 °), C22 = 9.8 x 10- 3 e-^ 149 °), (M° u , M° d ) « (100, 1.1) GeV. 

For the sake of simplicity and economy, we have set Q 2 = in this fit; however, values of 
I C22 1 ~ 10~ 3 can l ea d to similar results. Note that the magnitudes of the real parts of r], cr, e, 
and e' are nearly the same as those given in the CP-conserving case [9] (see Eq. (JBJ)); in 
particular the relative signs of these real parts are identical. The fit shown above leads to 
the following values for the fermion masses and mixings, while preserving the predictions for 
the neutrino system as in Eq. (JJj): 

(mf ys , m b (m b ), m T ) w (174, 4.97, 1.78) GeV 
(m c (m c ), m s {lGeV), m M ) « (1320, 101, 109) MeV 

(m° u (lGeV), m° d (lGeV), m°) « (10.1, 3.7, 0.13) MeV (21) 
{V us , V cb , \V ub \, \V td \)(<m z ) « (0.2250, 0.0412, 0.0037, 0.0086) 
= 0.150, 77^ = 0.374 

In obtaining the fermion masses at the low scales, we have not directly used Mf^ and 
M. d of Eq. (j3J). Instead, we have used: (a) m t (m t ) = 167 GeV and m T (m T ) = 1.777 GeK 
as inputs; (b) the GUT-scale predictions of our model for the ratios of masses - such as 
m b /m T , m UyC /m t , m^ s jm b ^ m^/m T etc; (c) renormalization in 2-loop QCD of these ratios in 
going from the GUT-scale to an effective SUSY-scale Ms = 500 GeV; and (d) the evolutions 
in 3-loop QCD and 1-loop QED of individual fermion masses as the effective momentum 
runs from Ms to the appropriate low energy scales [27] 12 . 

The primes on p' w and fj' w signify that these values are obtained from the SO(10)/G(224) 
model based fermion mass matrices (as in Eq. (j3j)), in conjunction with fermion masses and 
neutrino oscillations, as opposed to SM-based phenomeno logical values (Eq. (J2J)). 

12 Defining n a / b = (m a /m b ) GUT /(m a /m b ) Ms and r\ s = ■m f (M s )/m f (^ low ), wc get (for tan/3 = 5 and 03 (M z )=0.118): 
r] b/T = 0.6430, r) UjC/t = 0.4456, r; d>s/i , = 0.7660, T] et)l/T = 0.9999, r] u = 0.3954, rj d , s = 0.3982, Vc = 0.4418, rj b = 
0.6053, ?j e , M = 0.9894, r/ T = 0.9914, r) t = 0.9427. The CKM elements at low scales are given by V a/3 (< m z ) = V a/3 (GUT)/ K a0 , 
where K a p fa 0.91 for a/3 = ub, cb, td, and ts and K a p 1 for the other elements. 



18 



Note that, except for the very light fermion masses (m°, m° d , and m°) which would need 
corrections of order 1 to few MeV [28], all the other quark-lepton masses and especially the 
CKM mixings are in good accord with observations (see values quoted below Eq. (J2J) or 
Ref. [29]), to within 10%. (As alluded to before, we should not of course expect the very 
light fermion masses to be described adequately by the gross pattern of the mass-matrices 
exhibited in Sec. 2. In particular the "11" entries in Eq. (j3J) (expected to be of order 
10~ 4 — 10~ 5 ) arising from higher dimensional operators, which have been neglected in Sec. 
2, can quite plausibly lead to a needed reduction in m u by about 6-8 MeV and an increase 
in (m e and rrid) 13 by nearly (0.36 and 2-3) MeV respectively, at the 1 GeV scale, without 
altering the CKM mixings). 

The important point is that the SO(10)/G(224)-model presented in Sec. 3 has turned out 
to be capable of yielding values for p' w and fj' w that are close to the SM values as desired, 
while simultaneously being able to yield fermion masses of the two heavy families, all the 
CKM elements and neutrino oscillations (see Sec. 2), in good accord with observations. This 
in itself is non-trivial. 

Before presenting the results for CP and flavor violations some comments are in order as 
regards the parameters of the model versus its predictiveness. As expected, introduction of 
(in general four) phases in the Dirac mass matrices clearly increase the number of parameters 
compared to that for the CP-conserving case [9]. As a result, as long as we confine to the 
realm consisting of (a) the fermion masses and mixings, (b) CP and flavor violations induced 
only by the SM interactions, and (c) neutrino oscillations, the predictiveness of the model 
is reduced considerably (compare with the CP-conserving case of Ref. [9], see Sec. 2), the 
number of parameters now being comparable to the number of observables. Nevertheless, 
some gross features of the predictions in fact survive, even in the realm mentioned above, 
simply because: (a) the entries in the mass-matrices, governed by flavor symmetries, are 
hierarchical with a pattern as in Eq. (jSJ); (b) the phases are constrained 14 to lie between 
to 27r, and, (c) the system itself is constrained by the group theory of SO(10)/SU(4) C . One 

13 Note that the "11" entry for the up sector can differ from that for the down sector even in sign because of contribution 
through the operator 16il6il6^(16£f/M")(S/M) n which contributes only to m d and m e (so that Sm d = Sm e at Mqut) 
but not to m u . 

14 For instance, consider the familiar relation \V U s\ = \y md/ms — e"^ \/m u /m c \ , that holds for a hierarchical pattern. Given 
^/m d /m s 0.22 and \/m u /m c as 0.07, we cannot of course predict V U s precisely without knowing the phase angle <f>. Yet, 
since <j> can vary only between to 2n, \V US \ must lie between 0.15 and 0.29 
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can argue that these features in turn pretty much ensure the gross nature of the following 
predictions: (i) m° b /m° T « 1, (ii) \V US \ ~ 0.2, (iii) \V ub \ « ^m u /m c |V c6 |, (iv) \V t d\ ~ 
y/m d /m 8 \V cb \, (v) m^/m^ ~ 1/10, (vi) ~ 1/10 eV, and, (vii) sin 2 29°^ « (0.8-0.99), 
despite large variations in the parameters. 

The real virtue of the model (including the phases) emerges, once one includes SUSY CP 
and flavor violations. In this case, the realm of observables and thereby the predictiveness of 
the model expands enormously. The set of observables now includes not only the four entities 
listed in Eq. ((TJ-i.e. , Am^, ex, ^ m B d and S(Bd —> J/ipKs)— but also a host of others, for 
which the predictions of the G(224)/SO(10) model including (SM'+ SUSY) contributions, 
can a priori differ significantly from those of the SM contributions. In particular, the set 
includes observables such as (v) e' K (vi) Am^, (vii) S(B d — > 4>Ks), (viii) S(B d — > rj'Ks), 
(ix) S(B S — ► J/ip<f>), (x) S(B S — > (pKs), (xi) B — > Kir, (xii) B — > im (rates and asymmetry 
parameters), (xiii) b — > sj, (xiv) electric dipole moments of (n,e,Hg,d) and (xv) Lepton 
flavor violating processes (/i — > cy, r — > fij, r — > ej), and more. 

Now, the SUSY contributions do of course depend in part on the flavor preserving SUSY- 
parameters (i.e. m , m^, /i, and tan/5; we set A Q = at M*). But these few parameters 
should be regarded as extraneous to the present model, and hopefully, they would be deter- 
mined through the discovery of SUSY at the LHC The interesting point is that for a given 
choice of these flavor-preserving SUSY parameters (essentially m and m 1( / 2 ) the SUSY con- 
tributions to all the CP and/or flavor- violating processes listed above get completely deter- 
mined within our model, in magnitude as well as in phases. This is because all the flavor 
and in general CP violating sfermion (mass ) 2 -parameters 15 {{8m 2 ) l [ L RR LR ), arising through 
SO(10)/G(224)-based RG running from M* to Mqut are completely fixed in the model in 
terms of the parameters of the fermion mass-matrices (see Eqs. fllHj) . (fT%|) and (j!9)l ). The 
latter are, however essentially fixed by fermion masses and mixings, as shown in the fit given 
above (Eq. (j2~Uj) ). especially when we demand that the p' w and fj' w be close to the SM- values. 
In short, the inclusion of SUSY CP and flavor violations, treated in conjunction with fermion 
masses and neutrino oscillations, encompasses a host of processes without introducing new 
parameters and thereby increases the predictiveness of the model enormously. 

In this paper, we will present the results for some of the processes listed above, in par- 

15 counting the number of such (mass) 2 -parameters 
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ticular those shown in Eq. (JTJ) as well as those for e' K and the edm's of the neutron and the 
electron. Some of the other processes including lepton flavor violation will be considered in 
a separate paper. 

Using Eqs. (fT^j) and (fT3j) - (fT^j) for the squarks (mass) 2 elements (^Ilrrlr e ^ c as 
predicted in our model, the expressions given in Refs. [30-33] for the SUSY contributions, 
and the values of m s , m c and the CKM elements (including p' w and fj' w ) as obtained in the fit 
given above (see Eqs. (f2T)|) and (J2IJ)), we can now derive the values of the four entities listed 
in Eq. (JTJ), treating separately the cases of the SO(10) and the G(224)-models. For reasons 
explained below Eq. (JUJ), the SUSY contributions are reduced (in most cases) by about 
a factor of 2.5 in the amplitude for the case of G(224) compared to that of SO(10), being 
the effective symmetry in 4D. This distinction, as we will see, provides a way to distinguish 
between the SO(10) and the G(224)-models experimentally. The predictions of the model 
(corresponding to the fit shown in Eq. ()2()j) ) are shown in table 1. We have included both the 
SM' and the SUSY contributions in obtaining the total contributions (denoted by "Tot"). 
In quoting the numbers we have fixed, for concreteness, M* /Mqut ~ 3 and thus £ m 0.4, 
and have made a plausible choice for the SUSY spectrum - i.e. m sq « (0.8 - 1) TeV with 
x = m 2 g/m 2 sq ~ 0.8, although a variation in these parameters with m sq as low as about 600 
GeV or x = 0.5 — 0.6 can still lead to the desired results for all four quantities especially for 
the case of G(224) (see remarks below). 



(m , m 1/2 )(GeV) 


(800, 250) 


(600, 300) 




(a) 
SO(10) 


(b) 
G(224) 


(c) 
SO(10) 


(d) 
G(224) 


Am^-(Tot « SAT') (GeV) 


2.9 x 10~ 15 


2.9xl0- 15 


2.9xl0~ 15 


2.9X10" 15 


e K (SM') 


2.83xl0~ 3 


2.83x10-3 


2.83xl0- 3 


2.83xl0~ 3 


e K {Tot) 


1.30xl0 -3 


2.32xl0- 3 


2.01 xl0~ 3 


2.56xl0- 3 


Am Bd (Tot » SAT) (GeV) 


3.62xl0~ 13 


3.56xl0~ 13 


3.58xl0- 13 


3.55xl0- 13 


S(B d -> J/ipRs) (Tot « SM') 


0.740 


0.728 


0.732 


0.726 



Table 1. Predictions of the SUSY SO(10) and G(224) models corresponding to the fit for the fermion mass- 
parameters shown in Eq. (I20H . Either model with the fit as in Eq. 121)1 leads to the fermion masses and CKM mixings 
in good agreement with the data (see Eq. 1211 1. The total contribution (denoted by "Tot") represents the sum of 
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the SM' and the SUSY contributions. Note that the SUSY contribution is important only for ejf, furthermore it is 
relatively negative (as desired) compared to the SM' contribution (see discussion in text). The superscript s.d. on 
ArriK represents short distance contribution. 

In obtaining the entries for the i^-system we have used central values of the matrix element 
Bk and the loop functions T]i (see Refs. [1,33] for definitions and values) characterizing short 
distance QCD effects - i.e. B K = 0.86 ± 0.13, f K = 159 MeV, r) x = 1.38 ± 0.20, 16 rj 2 = 
0.57±0.01 and 773 = 0.47±0.04. For the -B-system we use the central values of the unquenched 
lattice results: / Bd v/^ = 215(11)(^ 3 )(15) MeV and fs s \Jb~~b s = 245(10) (^)(+£) MeV 
[34] . Note that the uncertainties in some of these hadronic parameters are in the range of 
15%; thus the predictions of our model as well as that of the SM would be uncertain at 
present to the same extent. 

At this stage the following comments are in order. 

(1) In the cases of Am^, Am^ and S(Bd — > J/tpKs), the SUSY contributions (with 
m sq ~ 0.8-1 TeV and x ~ 0.5-0.8) are found to be rather small (~ 0.5%, 2%, and 3% 
respectively) compared to the SM' contribution. As a result, for these three entities, the 
SM' contribution practically coincides with the total contribution, which is what is shown 
in the table. By contrast, for the same spectrum, the SUSY- contribution to ex is found to 
be rather sizeable (~ 20-25% ) 17 , and importantly enough, negative compared to the SM'- 
contribution 18 . The fact that it is relatively negative is an outcome of the model and, as it 
turns out, is most desirable (see below). 

(2) Comparing the predicted values shown in Table 1 with the observed ones (see those 
listed below Eq. ©), together with Am* = 3.47 x 10~ 15 GeV, we see that all four entities 
including ex and the asymmetry parameter S(B d — > J/ipK s ) agree with the data quite well, 
for the cases of SO(10) as well as G(224) shown in the last two columns (i.e. for m sq ~ 1 
TeV, and x m 0.8), and also for the case of G(224) in the second column (m sq m 800 GeV, 
x ~ 0.8). In making this comparison we are allowing for plausible (at present theoretically 

16 We will be guided by the error of ±0.20 on rji, used in [33], although that quoted in [1] is considerably larger (± 0.53). 
17 The fact that the SUSY contribution to ex (in contrast to those for Amj-, Ams d and S(B,i — » J/ipKg)) is relatively large 

is simply because the SM contribution to eg is strongly suppressed owing to the smallness of the relevant CKM mixings. 
18 In as much as we require p' w and fj' w to be close to the SM-based phenomenological values (as in Eq. 0), in accord with 

the observed values of the fermion masses, CKM-elements and neutrino oscillation parameters, we find that the class of fits 

satisfying these requirements lead to SUSY-contribution to ex that is relatively negative compared to the SM'-contribution. 
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uncertain but allowed) long distance contribution to Am^(~ ±15%), and uncertainties in 



B K or T]x < 10% (see entries for e K in the last three columns) and that in fB d \J Bs d by about 



(3) We note that a choice of the SUSY-parameters, e.g. (m , mi/ 2 ) = (800, 250) GeV, 
shown in the table, is in accord with the WMAP-constraint on CMSSM-spectrum in the 
event that the lightest neutralino is the LSP and represents cold dark matter [35]. 

(4) It is crucial that the SUSY contribution to ex (as mentioned above) is significant and is 
negative relative to the SW-contribution. Indeed this is what makes it possible for ex {Tot) 
to be desirably lower than the 6k(SM') = 2.83 x 10~ 3 and thereby to agree better in the last 
three columns (cases b, c and d) with e°^ s = 2.27 x 10~ 3 . Had the SUSY contribution been 
positive relative to the SM' contribution and still as significant as above, ex {Tot) would 
have been (3.34, 3.53, and 3.10) xlO -3 for the cases (b), (c), and (d) respectively, in strong 
disagreement with observation. In short, the SUSY contribution of the model to €k has just 
the right sign and nearly the right magnitude to play the desired role. This seems to be an 
intriguing feature of the model. 

(5) Since the values of the CKM elements including p' w and ff w obtained within our model 
(see Eq. (|2*T|) ) are quite close to the SM based phenomenological values (see Ref. [1] and 
Eq. 0), we would of course expect that the SM' contributions should nearly be the same 
as the SM contributions, for the same choice of the hadronic parameters {Br, T)i, fs d \J BB d 
etc.). For instance, using the central values of the parameters given in the recent update 
of the CKM-triangle analysis by M. Bona et al. [2], that is, A = \V US \ = 0.2265, \V c b\ = 
4.14 x 10" 2 , p w = 0.172, fj w = 0.348, m c = 1.3 GeV and f K = 159 MeV, and the hadronic 
parameters as in our case - that is, B K = 0.86, 771 = 1.38, r] 2 = 0.57, and r/ 3 = 0.47 - one 
obtains €k{SM) = 2.72 x 10~ 3 which is about 20% higher than the observed value. Contrast 
this with the predictions for (Tot) of the SO(10) or G(224) models for the cases (b), (c) 
and (d) in Table 1 where the discrepancies between the predicted and observed value of ex 
range from 2 to 12% with varying signs. At present, such discrepancies, even as high as 20% 
for the SM, can of course be accommodated by allowing for uncertainties in B K , r]i, and 
also in A. 

(6) One main point we wish to stress here, however, is this: At present, the distinctions 
between the predictions of the SM (in particular for ex) versus those of the SUSY SO(10) 
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or G(224) models on the one hand, and those between the predictions of the SUSY SO(10) 
versus the G(224) models on the other hand (compare columns (a), (b), (c) and (d) of Table 
1) are marred in part because of uncertainties (~ 15%) in the hadronic parameters (Bk, f]i 
etc.) as well as that (~ 2%) in A, and in part because SUSY is not discovered as yet, and thus 
the SUSY spectrum is unknown. But once (hopefully) SUSY is discovered at the LHC and 
thereby the SUSY parameters get fixed, and in addition the uncertainties in the hadronic 
parameters are reduced (hopefully) to a few percent level through improvements in the lattice 
calculations, we see from the analysis presented above that we can utilize the combined set of 
informations to distinguish experimentally between the SM versus the SUSY SO(10)-model 
versus the SUSY G(224)-model. It is intriguing to see that even low energy experiments 
involving CP and flavor violations can help distinguish between the SO(10) versus the G(224) 
models, both of which nearly coincide as regards their predictions for fermion masses and 
neutrino oscillations. In this way they can shed light on physics at the superheavy scale 
M* > Mqut- The experimental distinctions will of course be even sharper once we include 
predictions for the other processes, some of which are presented below. 

(7) Bj — > <pKg, AmB s and b — ► 57: We now consider the CP violating asymmetry 
parameter S(Bd — > 4>Ks). For a representative choice of (m , mi/2) = (600, 300) GeV, we get 
5f L = (1.40-0.012i)xl0- 2 , 5f R = -(5.39+6.27i)xl0" 3 , 5f R = -(0.29+3.08i) x 10~ 4 /tan/3 
and 5 2 R \ = -(1.92 + 2. 70i) x 10~ 4 /tan/3 as predictions of our model (see Eqs. (JI3J) and (fP3])). 
It is easy to verify that the SUSY-amplitude for this decay in our model is only of order 1% 
(or less) [36] compared to that in the SM. As a result, adding SM' and SUSY contributions 
to the decay amplitudes, we obtain: 

S(B d ^ 4>K S ) (Tot » SM') » 0.728 (22) 

Allowing for variant fits which also give fermion masses and CKM mixings in good agree- 
ment with observations, we find that S(Bd —> <pKs) should lie in the range of ~ +0.65 
to +0.73. Thus our model predicts that S(Bd — > (j)Ks) is close to the SM prediction 
(ph 0.70 ± 0.10) and certainly not negative in sign. 19 When we started writing this pa- 
per, BaBar and BELLE data were yielding widely differing values of (0.45 + 0.43 + 0.07) and 
(—0.96+0. 50l ; 07 ) respectively for S(Bd — > 4>Ks) [37]. Most recently, the two groups reported 

19 Our prediction in this regard was reported at the Fujihara seminar [7], held February 23-25, 2004. 
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new values for the asymmetry parameter S(E>d — > 4>Ks) = (+0.50±0.25toM)BaBar', (+0.06 ± 
0.33 ± 0.09) belle [37], at the Beijing International Conference on High Energy Physics . 
Meanwhile there have been many theoretical and phenomenological attempts [36,38] to ob- 
tain possible large deviations in S(Bd — > <pKs) from the SM-value, including, in some cases, 
negative values for the same (as suggested by the earlier BELLE data). It will thus be ex- 
tremely interesting from the viewpoint of the G(224)/SO(10)-framework presented here to 
see whether the true value of S(B d — > 4>K S ) will turn out to be close to the SM-prediction 
or not. 

Including contributions from 8 2 L \ RR and 5 2 L \ RL ( as predicted in our model), we get: 



f \ B 

Am Bs (Tot « SM') » 17.3 P^i ^^ T ■ (23) 

This is of course compatible with the present lower limit on Am# s > 14.4p.s~ 1 [2]. 

Using 5 2 RL given above, we obtain A{pL — > s R j)g ~ (1 — 1.5) x 10~ 10 GeV^ 1 / tan/3. Even 
allowing for variant fits, the SUSY-amplitude, in this case, is found to be only about (1.5- 
5)% of the SM amplitude. The same conclusion holds also for A(1)r — > SLj)g- In short, our 
results for (B d — > (f>Ks), Am^ s and b — > nearly coincide with those of the SM. 

(8) Contribution of the A term to e' K : Direct CP violation in Kl — > vr7r receives a new 
contribution from the chromomagnetic operator Q~ = (g /lQT[ 2 ){si J a^ u t a d R — s R a^ u t a di J )G a llu) 
which is induced by the gluino penguin diagram. This contribution is proportional to X 2 \ = 
Im[{5 ( [ R )2i — {&lr)i2\-> which is predictable in our model (see Eqs. (|TR|) and (fTHjl ). Following 
Refs. [39] and [40], one obtains: 

n, o 410 Mey w 500 GeF w fnA . 

Re(e'/e) § w 91 B G ( ■ ) ( )X 21 24 

m s + ma m~ g 

where Bq is the relevant hadronic matrix element. Model-dependent considerations (allowing 
for m 2 K /m 2 corrections) indicate that Bq ~ 1-4, and that it is positive [39]. Using the 
prediction of our model (via Eqs. (fTHj) and (fT9*j)). for a typical SUSY- spectrum used in 
previous considerations (e.g. (m G , mi/2) = (600, 300) GeV), we obtain: X 21 w 2.1 x 
10~ 5 / tan (3. Note that the sign of X21, as derived in the model, is positive. Inserting this in 
Eq. (J2U), and putting (m s + m d ) m 110 MeV, we get: 

Re(e'/e)~ g w +(8.8 x 10- 4 )(S G /4)(5/tan^) . (25) 
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We see that if the positive sign of Bq is confirmed by reliable lattice calculations, the gluino 
penguin contribution in our model can quite plausibly give a significant positive contribution 
to Re(e'/e)g « (4—14) x 10~ 4 , depending upon B G m 2—4 and tan/? m (3—10). At present the 
status of SM contribution to Re(e' / e) is rather uncertain. For instance, the results of Ref. [41] 
and [42] based on quenched lattice calculations in the lowest order chiral perturbation theory 
suggest negative central values for Re(e'/e). (To be specific Ref. [41] yields Re(e'/e)sM = 
(—4.0 ± 2.3) x 10~ 4 , the errors being statistical only.) On the other hand, using methods 
of partial quenching [43] and staggered fermions, positive values of Re(e'/e) in the range of 
about (3 — 13) x 10 -4 are obtained in [44]. In addition, a recent non-lattice calculation based 
on next-to-leading order chiral perturbation theory yields Re(e'/e)sM — (19±2lg±6) x 10~ 4 
[45]. The systematic errors in these calculations are at present hard to estimate. The point 
we wish to note here is that the SUSY-contribution to Re(e'/e), in our model, is significant, 
and when the dust settles, following a reliable calculation of Re(e'/e) in the SM, it would be 
extremely interesting to check whether the SUSY-contribution obtained here is playing an 
important role in accounting for the observed value given by Re(e'/e) b s = (17±2) x 10~ 4 [46] 
or not. 

(9) EDM of the neutron and the electron: RG-induced A-terms of the model gen- 
erate chirality-flipping sfermion mixing terms (S£^)ij, whose magnitudes and phases are 
predictable in the model (see Eqs. (fT£j) and (JE3)), for a given choice of the universal SUSY- 
parameters (m , mu2, and tan/?). These contribute to the EDM's of the quarks and the 
electron by utilizing dominantly the gluino and the neutralino loops respectively. We will ap- 
proximate the latter by using the bino-loop. These contributions are given by (see e.g. [47]): 

(d d ,d u ) Aind = (-11) f e ^/(^/mftn 

m - rn 2 - I 26 ) 

For a representative choice (m ,mi/ 2 ) = (600, 300) GeV (i.e. m sq = 1 TeV, m- g = 900 GeV, 
mj = 636 GeV and m s = 120 GeV ), using Eqs. (fTS jl and (fTH jl . we get: 



4.15 x 10~ 26 ,, N , _ _ ofi , 1.1 x 10" 28 



ecm 



(ddiA ^ = — 7. ecm: ((L)a' j = (—1-6 x 10 ) ecm: (d e ) A . , 

v ' md tan/3 v 1 md y 1 v ' md tan/3 

The EDM of the neutron is given by d n = \(Add — d u ). Thus for SO(10), with the choice of 

(m D , mi/2) as above, we get 

(d n ) Amd = (1.6, 1.08) x 10" 26 ecm for tan/? = (5, 10) . (28) 
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Note that these induced A-term contributions are larger for smaller tan/3. For an alter- 
native choice (m , mi/2) = (800, 250) GeV, which as mentioned before is compatible with 
the WMAP/CDM-constraint [35], the predicted EDM of the neutron and the electron are 
reduced respectively by about 36% and a factor of 3.6. The predictions for the G(224)-model 
are smaller than those for the SO(10)-model by nearly a factor of two in all cases. 

We should also note that intrinsic SUSY phases denoted by (ft a = o J rg(Alm 1 / 2 ) and 
4>b — < ir g{ m i/2^B*), if present, would make additional contributions to EDM's through 
gluino and/or chargino/neutralino loops, which should be added to the contributions shown 
above. These contributions have been widely discussed in the literature (see e.g. [47]). As is 
well known, with A a = or small (< 1 GeV) at the scale M* (as we have chosen, following the 
examples of Refs. [19] and [21]), these contributions are proportional to (m^e)/^ tan /3(sin </>#). 
They would be typically about 50-300 times larger than the values shown above (Eqs. (J27|) 
and (J28|) ). if the relevant intrinsic SUSY phases are nearly unity. This is the familiar SUSY 
CP problem. The point of the present paper is that even if the intrinsic SUSY phases are 
naturally zero or insignificantly small, as would be true in a theory where the SUSY CP 
problem is naturally solved 20 ' 21 , the induced A-term contributions to EDM's (arising from 
GUT-scale physics) as presented above would still be present. The interesting point is that 
these contributions are completely determined in magnitude and phase within our model (for 
a given choice of the SUSY universal parameters (m G , m^, tan/?)). 

Given the experimental limits d n < 6.3 x 10 -26 e cm [51] and d e < 4.3 x 10 -27 e cm [52], we 

20 A possible solution to the SUSY CP problem could arise as follows. Assume that CP violation arises spontaneously, only 
in the visible sector, through the VEV's of fields at the GUT-scale, like those of 16h, 45 h and the singlet S. One can 
argue that the VEV's of at least some of these can be naturally complex or purely imaginary, consistent with the minimization 
of the potential, even if all parameters in the potential are real. In this case, intrinsic SUSY phases like those in m^2, A and B 
are, of course, zero. Now if the /j-term can be derived, through a satisfactory resolution of the /^-problem, for example, either 
by the Giudice-Masiero mechanism [48], or by the ideas suggested in [49], or by involving a coupling in the supcrpotcntial of 
the form [50]: KlO/jlO//A r + XN [i + where the singlet N is not allowed to couple to the other fields mentioned above, and 
acquires a real VEV of order 1 TeV (as needed), with k and A being real, then the ^t-term would also be real. In this case, all 

intrinsic SUSY phases would disappear. We plan to explore this possibility in a future work. 

21 An alternative resolution of the SUSY CP problem arises in a class of gaugino mediated SUSY-breaking (with the fi- 

problem being resolved for example as in [48]) in which all relevant SUSY parameters become proportional to m 1( i 2 [19]. A 
third resolution of the SUSY CP problem would arise in a model where both A and B terms are naturally zero or sufficiently 
small at the scale M*. This is precisely what happens in, for example, the anomalous U(1)a model of SUSY-breaking that 
arises in the context of a three-family string-solution [21]. In this case, extra gauge symmetries of the model suppress both A 
and B terms at M* . 
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see that the predictions of the model (arising only from the induced A-term contributions) 
especially for the EDM of the neutron is in an extremely interesting range suggesting that 
it should be discovered with an improvement of the current limit by a factor of about 10. 

6 Summary and Concluding Remarks 

In this paper we have explored the possibility that (a) fermion masses, (b) neutrino oscilla- 
tions, (c) CP-non conservation and (d) flavor violations get intimately linked to each other 
within supersymmetric grand unification, especially that based on the symmetry SO(10) or 
an effective symmetry G(224)= SU(2) L x SU(2) R x SU(A) C . In this context, we extend 
the framework proposed previously in [9], which successfully described fermion masses and 
neutrino oscillations (see Sec. 2), to include CP violation. We assume, in the interest of pre- 
dictiveness, that CP-violation, arising through the SM as well as SUSY interactions, has its 
origin entirely (or primarily) through phases in the fermion mass matrices. We also assume 
that flavor-blind universal SUSY parameters (m D , and tan/3 with A Q being small or real) 
characterize SUSY-breaking effects at a high scale M* > M GUT . In this case, all the weak 
scale CP and/or flavor-violating as well as flavor-preserving sfermion transition- elements 
^llrrlr' an d the induced A-parameters, get fully determined within the model, in their 
magnitudes as well as in phases, simply by the entries in the SO(10)-based fermion mass- 
matrices, once the few soft parameters (m G , and tan/3) are specified. This is how CP 
and flavor violations arising jointly from the SM and SUSY interactions, get intimately tied 
to fermion masses and neutrino oscillations, within a predictive SO(10)/G(224)-framework 
outlined above. The presence of GUT-scale physics induces enhanced flavor violation with 
and without CP violation, which provides a distinguishing feature of the model 22 , relative to 
CMSSM or MSSM. 23 

22 Even in the case of CMSSM, all the parameters of MSSM at the electroweak scale (some 105 of them) are of course also all 
fully determined in terms of the SUSY-parameters (m , m 1 / 2 and tan/3) and those involving the fermion masses and mixings. 
However, in this case, as mentioned in Sec. 4, owing to the absence of GUT-scale physics in the interval M* — » Mqut, the 
most interesting effects on the entities considered here (e.g. those on ek, e' K an d tne EDM's) would be absent or negligibly 
small. 

23 While we have focusscd in this paper on the SO(10)/G(224)-model of Ref. [9], we note that generically such enhanced 
flavor and/or CP violations arising from GUT-scale physics would of course be present in alternative models of SUSY grand 
unification [53] as well, as long as the messenger scale M* lies above Mqjjt- The detailed predictions and consistency of any 
such model as regards flavor and/or CP violations can however depend (even sensitively) upon the model, and this is worth 
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As mentioned in Sees. 1 and 5, the framework presented above faces, however, a prima- 
facie challenge. Including SM and SUSY contributions, the question arises, can the frame- 
work successfully describe the observed features of CP and flavor-violations including those 
listed in Eq. (0), while retaining the successes of the CP-preserving framework [9] as re- 
gards fermion masses and neutrino oscillations? Our work here shows that the SUSY 
SO(10)/G(224)-framework proposed here, which extends the framework of Ref. [9], indeed 
meets this challenge squarely. In the process, it makes several predictions, only some of 
which are considered here; these can eventually help distinguish the framework from other 
alternatives. 

Our results can be summarized as follows: 

(1) It is found that, with allowance for phases, there exists a fit to the parameters of the 
fermion mass-matrices (Eq. (j^J) which successfully describes fermion masses, all the CKM 
elements and neutrino oscillations as in Ref. [9] (see Eq. (|21Jl ). and simultaneously yields the 
Wolfenstein parameters (p' w , fj' w ) that are close to the phenomenological SM values (Eq. 
(0)). The merit of obtaining such values for (p' w , fj' w ) in accounting for the data on CP and 
flavor violations in quantities such as those listed in Eq. ([TJ has been stressed in Sec. 5. 

(2) With these values of (p' w , fj' w ), and a plausible choice of the SUSY-spectrum 24 (i.e. 
m sq ^(600-1000) GeV and rrig ^(500-900) GeV, say, it is found that the derived values of 
all four quantities (i.e. Arrix, £k, A.rriB d and S(Bd — > J/ipKs) agree with the data quite 
well (allowing for up to 15% uncertainty in hadronic matrix elements), see Table 1. 

(3) Although the SM' contributions (for the fit shown in Eq. ([20))) nearly coincide with 
the SM contributions to all entities, and SUSY-contributions to entities such as Am#, Am# d 
and S(Bd — > J/ipKg) are rather small (^ a few%), the contributions from SUSY, as a rule, 
are nevertheless prominent (of order 20-25%) especially when the SM (or SM') contributions 
are suppressed (for example due to smallness of the mixing angles). Such is precisely the 
case for sr, e' K and the edm's of the neutron and the electron, (as well as for lepton flavor 
violating processes to be discussed in a forthcoming paper [25]). It is found that the SUSY- 
contribution to is sizable (of order 20-25%) and negative relative to the SM' contribution, 
just as desired, to yield better agreement between the predicted and the observed value (see 

examining. 

24 Lighter masses for the SUSY particles like m sq ?s 600 GeV and rrig PS 500 GeV (say) are allowed for the case of G(224), 
though not for SO(10) (see discussion following Table 1). 
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Table 1). 

(4) The sizable negative contribution of SUSY to ek in our model provides an important 
tool to help distinguish not only between the SM versus the SUSY SO(10)/G(224)-models, 
but also between the SO(10) and the G(224)-models themselves (see Table 1). Such distinc- 
tions would be possible once (hopefully) SUSY is discovered at the LHC and thereby the 
SUSY parameters get fixed, and in addition the uncertainties in the hadronic parameters 
(B K and 771) are reduced to (say) a 5% level or better, through improved lattice calculations. 

(5) The model predicts that S(B d — > 4>K S ) should lie in the range of +(0.65-0.73), i.e. 
close to the SM-prediction. Given the present still significant disparities between the BaBar 
and BELLE results versus the SM-predictions, it would be interesting to see where the true 
value of S(Bd — > 4>Ks) would turn out to lie. 

(6) It is interesting that the quantity X 2 \ = Im[( y 5'[ R )2i — {.^lr)i2\i relevant for e' K , is 
found to be positive in the model. If the presently indicated sign of the relevant hadronic 
matrix element Bq being positive is confirmed, our model would give a positive contribution 
to Re(e'/e) which quite plausibly can lie in the range of +(4 — 14) x 10~ 4 . While this is in 
the interesting range, its relevance can be assessed only after the associated matrix elements 
are determined reliably. 

(7) The model predicts that the EDM's of the neutron and the electron should be discov- 
ered with improvements in the current limits by factors of 10 and 100 respectively. (Intrinsic 
SUSY-phases, even if present, would not alter this conclusion as long as there is no large 
cancellation between different contributions.) 

(8) It would be most interesting to explore the consequences of the model, involving 
SUSY contributions, to other processes such as B s — > J/ip(fi, B s — > (pKs, B — > Kn, B — > nn, 
B — > DK, K L — > ttuu, K + — > nuu, and especially lepton violating processes (such as 
fj, — > e7, t — > /ry, r — > cy etc.). We stress that the net (SM' + SUSY)-contributions to all 
these processes involving CP and/or flavor violations are completely determined within our 
model. They do not involve any new parameters. For this reason, the model turns out to 
be highly predictive and thoroughly testable. These processes are under study; a paper on 
lepton flavor violation is in preparation. 

To conclude, the SUSY SO(10)/G(224) framework, as proposed in Ref. [9] and extended 
here, subject to the assumption of universality of SUSY parameters, drastically reduces the 
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parameters for SUSY-contributions to CP and flavor- violations. In effect, the extension 
proposed here ties in fermion masses, neutrino oscillations, CP and flavor violations within 
a predictive and testable framework. 
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